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Phonon coupling (PC) corrections to magnetic moments of odd neighbors of magic and semi- 
magic nuclei are analyzed within the self-consistent Theory of Finite Fermi Systems (TFFS) based 
on the Energy Density Functional by Fayans et al. The perturbation theory in g\ is used where 
pi, is the phonon-particle coupling vertex. A model is developed with separating non-regular PC 
contributions, the rest is supposed to be regular and included into the standard TFFS parameters. 
An ansatz is proposed to take into account the so-called tadpole term which ensures the total 
angular momentum conservation with g\ accuracy. An approximate method is suggested to take 
into account higher order terms in g\. Calculations are carried out for four odd-proton chains, the 
odd Tl, Bi, In and Sb ones. Different PC corrections strongly cancel each other. In the result, the 
total PC correction to the magnetic moment in magic nuclei is, as a rule, negligible. In non-magic 
nuclei considered it is noticeable and, with only one exception, negative. On average it is of the 
order of —(0.1 -f- 0.5) /j,n and improves the agreement of the theory with the data. Simultaneously 
we calculated the gyromagnetic ratios </£ h of all low- lying phonons in 208 Pb. For the 3]~ state it is 
rather close to the Bohr-Mottelson model prediction whereas for other L-phonons, two 5~ and six 
positive parity states, the difference from the Bohr-Mottelson values is significant. 

PACS numbers: 21.60.Jz, 21.10.Ky, 21.10.Ft, 21.10.Re 



I. INTRODUCTION 

This work is devoted to studying effects of the in- 
teraction of low- lying collective states ( "phonons" ) with 
single-particle degrees of freedom. The first microscopic 
description of phonons in spherical nuclei was made more 
than 50 years ago by Spartak Belyaev in the seminal ar- 
tide i. 

In the last decade, a considerable success of different 
approaches based on the mean field theory in the descrip- 
tion of nuclear bulk properties was achieved. It concerns 
the Hartree-Fock (HF) method with Skyrme force 2, 3], 
the HF method with Gogny force [J], the Relativistic 
Mean-Field approach Q and the so-called Generalized 
Energy Density Functional (EDF) method by Fayans et 
al. Q- The last approach is essentially close to the self- 
consistent Theory of Finite Fermi Systems (TFFS) Q 
based on TFFS by Migdal [8] supplemented with the 
Fayans-Khodel self-consistency relation 0. A specific 
feature of the Fayans method is the consideration of the 
pairing problem by solving the Gor'kov equations in the 
coordinate space with the method developed by Belyaev 
et al. [lC| . Recently, all these approaches were also suc- 
cessfully applied to the description of phonons in spher- 
ical nuclei within the self-consistent Quasiparticle Ran- 
dom Phase Approximation (QRPA) or similar methods, 



see e.g. ^MM- 

The next step consists in taking into account the 
phonon coupling (PC) effects. These effects were studied 
in detail in the Quasiparticle-Phonon Model by Soloviev 
[lH and, on the phenomenological level, within the so- 
called Nuclear Field Theory by Bortignon, Broglia et 
al. (l7j . On a more microscopic level, they were con- 
sidered mainly for describing different kinds of Giant 
Resonances and the Pygmy-dipole Resonance within the 
(Q)RPA+PC method and Extended Theory of Fi- 
nite Fermi Systems (ETFFS) [Hj]. Self-consistent ex- 
tensions of these methods have been proposed recently, 
for example, (Q)RPA+PC [20] and the ETFFS in the 
Quasiparticle Time Blocking Approximation (QTBA) 
"111 [ETFFS(QTBA)] [22j, see also the recent review 
23| devoted to the Pygmy-dipole Resonance. It should 
be mentioned as well successful applications of the self- 
consistent ETFFS (QTBA) in [HI and in [H[. The ap- 
proach is limited to magic and semi-magic nuclei where 
there is a small parameter , the square of the phonon 
creation amplitude. 

Here we concentrate on the self-consistent description 
of the PC corrections to magnetic moments. Due to mod- 
ern Radioactive Ion Beam facilities a lot of new data on 
magnetic and quadrupole moments appeared recently in- 
cluding nuclei distant from the beta-decay valley. The 
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bulk of the data obtained till 2005 is collected in a very 
comprehensive compilation by Stone [26j | . More recent 
data are presented in original articles, e.g. in Ref. [13]; 
where new, together with older, data on magnetic and 
quadrupole moments of a long chain of copper isotopes 
from N = 28 to N — 46 are presented and successfully 
analyzed within the Many-Particle Shell Model (MPSM) 
[28] | . This approach is very comprehensive as it takes into 
account all main inter-nuclear correlations. The necessity 
to introduce many parameters of the effective interaction, 
the single-particle mean field and the effective particle 
charges is some deficiency of the MPSM. In addition, the 
domain of the MPSM applications is, by technical rea- 
sons, limited to nuclei with A < 90 4- 100. 

For heavier nuclei, the challenge of experimentalists 
was partially responded within the self-consistent TFFS 
for mag netic moments [29l |30| and quadrupole moments 
Ull H2 of odd spherical nuclei. Mention also the 
first self-consistent calculation of the quadrupole mo- 
ments of the 2+ states [33[ . The calculations were mainly 
limited to semi-magic nuclei considered in the "single- 
quasiparticle approximation" where one quasiparticle in 
the fixed state A = (n, l,j, m) with the energy e\ is added 
to the even-even core. The QRPA-like TFFS equations 
for the nuclear response to the external field, magnetic 
in [29|, |30[ or quadrupole in [l5|, [3l|, HH , were solved on 
the base of the Energy Density Functional (EDF) by 
Fayans et al. @, [34| [35j]. For magnetic moments, the 
original Fayans functional DF3 @, |35| was used whereas 
for quadrupole moments it was used together with its 
modification, DF3-a, which was introduced in [36] to ex- 
tend this approach to nuclei heavier than lead. It differs 
from the original one by the spin-orbit and effective ten- 
sor terms which are important only for the fine structure 
of the single-particle spectrum. For the quadrupole mo- 
ments, the difference between the predictions of the two 
functionals turned out noticeable with a preference of 
DF3-a. We also mention that the same is true for the 
energies and excitation probabilites of the 2+ states in 
the lead and tin isotopes [TH, H3] ■ 

According to the TFFS [8] , a quasiparticle differs from 
a particle of the single-particle model in two respects. 
First, it possesses a local charge e q and, second, the core 
is polarized due to the interaction between the parti- 
cle and the core nucleons via the Landau-Migdal (LM) 
amplitude T . In other words, the quasiparticle pos- 
sesses the effective charge e c ff caused by the polarizability 
of the core which is found by solving the TFFS equa- 
tions. In the MPSM, a similar quantity is introduced 
as a phenomenological parameter. Thus, the consid- 
eration was made on the RPA (or QRPA) level which 
is the "zero approximation" to the problem, corrections 
due to phonon coupling (PC) effects were only estimated 
in [29M31I] . This does not concern the investigation of 
quadrupole moments of the 2\ states in Ref. [U as 
the effect under consideration there is beyond QRPA it- 
self. On average, reasonable description of the data was 
obtained in the Refs. cited above, with the accuracy 



of S[i ~ 0.1 -7- 0.2 n.m. for magnetic moments, and 
8Q ~ 0.1 -j- 0.2 b for quadrupole moments. This indi- 
cates that the single-quasiparticle approximation works 
well for such nuclei on average, and the PC effects are 
usually regular and included in the TFFS parameters, 
i.e. the charges e q and the LM amplitudes T . However, 
there are several "bad" cases, with &\x ~ 0.5 n.m. for 
magnetic moments, and 5Q ~ 0.5 b for quadrupole mo- 
ments which were attributed to some non-regular PC ef- 
fects. The estimations in [3(| for magnetic moments and 
in [3lj for quadrupole ones have shown that this inter- 
pretation looks reasonable and a more detailed analysis 
of the PC effects is necessary. 

In this context, it is worth citing older calculations of 
magnetic moments within the conventional TFFS with 
the use of the Saxon- Woods potential as the nuclear mean 
field MM- 



Dealing with the PC problem within TFFS or other 
microscopic approaches with fitted parameters, one faces 
the problem of refitting the parameters provided all the 
PC corrections are taken into account. Exactly this way 
was chosen in [24j |. where a new set parameters of the 



Skyrme functional was found which corrects for the PC 
effects. The same mode of action was used in the first ver- 
sion of the self-consistent TFFS 0- In this paper, dealing 
with magnetic moments of magic and semi-magic nuclei, 
we chose a simpler way by separating only the kinds of the 
PC diagrams which behave in non-regular way depending 
significantly on the nucleus and the single-particle state 
under consideration. The rest of the PC corrections, in 
agreement with the estimations, is supposed to be reg- 
ular and is included in the standard TFFS parameters. 
Within this model, we analyzed first odd neighbors of the 
magic 208 Pb nucleus. In this case the PC effects appear 
mainly due to the 3j~ state, although sometimes the sum 
of contributions of all other phonons is also important. 
For these nuclei, the problem was analyzed in the well- 
known paper by I. Hamamoto (4l| within the Nuclear 
Field Theory [l7|, the approach which operates with a 
set of phenomenological parameters for each nucleus un- 
der consideration. On a microscopic level, t his p roblem 
was considered by Platonov [42j and Tselyaev [43| . In the 
last article, some sets of higher order terms in g\ were 
summed up within QTBA. 

In non-magic nuclei, the contribution of the l\ state 
dominates. Again we limit ourselves to the semi-magic 
nuclei where this state is usually not too collective and 
the parameter g\ is, as a rule, small. In addition, we will 
consider only such odd semi-magic nuclei where the odd 
nucleon belongs to the normal sub-system which simpli- 
fies the calculations because all the equations for the PC 
corrections do not include pairing effects. For example, 
in the lead region we will analyze the odd Tl and Bi iso- 
topes, but not the odd Pb nuclei, except 207 > 209 pb. 
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FIG. 1: The effective magnetic field V(M1). The meaning 
of separate blocks becomes obvious after comparing with Eq. 

ra- 



il. SELF-CONSISTENT TFFS 

For completeness, we outline first the TFFS formalism 
for the effective magnetic field without phonon correc- 
tions. The magnetic moment of an odd nucleus or the 
probability of the Ml-transition in such a nucleus are 
determined from the effective field V which is related to 
the external field Vq = p,, where fi is the operator of the 
magnetic moment: 



(1) 



with gf=l, g?=0, 5 f=5.586, and 5 ™= - 3.826. The ef- 
fective field V obeys an equation which is similar to the 
RPA equation with the LM amplitude T playing the role 
of the effective NN-interaction. In systems without pair- 
ing, e.g. magic nuclei, this equation reads 

V(Ml;u) = Vo(Ml) + TA{lj)V{M\-uj), (2) 



where V (M1) 



e q (i and A{w) 



jG(s)G(e 



u>)ds / (2iri) stands for the particle-hole propagator, G(e) 
being the single-particle Green function. This is illus- 
trated diagrammatically on Fig. 1. The magnetic mo- 
ments of an odd nucleus with the odd nucleon in the 
one-quasiparticle state A = {n r ,l,j 1 m) — (v,m) are 
the diagonal matrix elements of the static effective field, 



fix = (A|V(Afl;w = 0)|A)|, 



whereas the matrix ele- 



ments of the Ml-transition Ai — !• A2 between two states 
with the energies £\ and e 2 respectively are given by 
the non-diagonal matrix elements M\ 2 = (Xi\V(Ml; lu = 
£2 - £i)|A 2 ). 

In the explicit form, the local quasiparticle charge e q 
in Eq. @ is defined as follows: 



I + (I - 20)r 3 



(9 p s + -<??)(! -2C s )r 3 



+ Ctl^n)©^] 1 . 



(3) 



Here, following Refs.[29|, |30(, in addition to the usual 
local charge parameters ( s ,(i of the FFS theory a 



new "tensor" (or "Z-forbidden" ) charge £t is introduced. 
Terms of similar structure appear naturally from MEC 
contributions to the magnetic moments 

In nuclei with pairing, Eq. ^ transforms to a set of 
three equations [8j 



V(u) = V r (w)+J 7 A(w)V(w), 



(4) 



where all the terms are matrices. In the standard TFFS 
notation Q, we have: 




Vq 




(5) 



where the fields dip denote variations of the gap func- 
tions AW't 2 ) in the external field Vq. The vector analogue 
of (O contains 3 <8> 3 matrices: 



and 



A(u) 



C{u) Mi{lo) M 2 (uj) 

o{lu) -MM N 2 {oj) 

0(-u) -M(-w) JV 2 (-w) 



(6) 



(7) 



where £, Aii, and so on stand for integrals over e of the 
products of different combinations of the Green function 
G(e) and two Gor'kov functios F^(e) and F^(e) 0. 

In the matrix F, the diagonal terms are the usual LM 
amplitude J- of the effective interaction in the particle- 
hole (ph) and the interaction amplitude J-^ irreducible in 
the pp (or hh) channel. The non-diagonal terms couple 
the ph channel with the pp and hh ones. As the analysis 
of [29( has shown, the interaction amplitude J-^ for the 
Ml symmetry is, evidently, small. Therefore, we put 
.7^=0. Then, it is natural also to put equal to zero all 
the non-diagonal terms of J 7 *. In this case, we have d\ = 
d 2 = and the only modification of Eq. ([2|) is the change 
A -> £, 



C(u)= / — G(e)G{e + u) - F 1 (e)F 2 (e + u)P 



de 

(8) 

where P = 1 for T-even fields and P = — 1 for the T- 
odd ones. To calculate the propagator C with consistent 
account for the continuum states, we use the generaliza- 
tion [45( for superfluid systems of the Shlomo-Bertsch 
method 

All calculations below will be carried out with the use 
of the self-consistent basis generated with the General- 
ized EDF by Fayans et al., 



Eq = J £ [p„(r),pp(r),j/„(r),i/ p (r)]cf 3 r, 



(9) 



depending simultaneously on the normal p r and anoma- 
lous v T densities. The DF3-a version of the normal EDF 



[36| will be used and the "surface" anomalous EDF Q. 
More details can be found in [15j . 

For magnetic moments we consider only the spin- 
dependent parameters contribute. The spin-dependent 
LM amplitude, just as in [29|, [3(J, is chosen in the form 
of 



^•spin j-spm 



(10) 



where the pion and /9-meson exchange terms are added 
to the central force term To- The central force was ap- 
proximated with the zero Landau harmonics. Following 
to [13, 13 , we take into account the momentum depen- 
dence of the amplitude g' ( "Migdal force" ) : 



-7-spin 



C [g (er lC r 2 ) + (<7i<7 2 )(tiT 2 )] , (11) 



/(?) = 



1 



(12) 



where the normalization factor Co — dn/deF) 1 = 
300 MeV- fm 3 , and the value of r = 0.4 fm was found in 
[47l . |48| for the parameter which determines the depen- 
dence of the amplitude g' on the momentum transfer. 
To evaluate the PC corrections we need in the vertex 
the creation amplitude of the L-phonon. It obeys 
the homogeneous equation corresponding to Eq. ([J), 



g L (u}) = FA{ui)g L {u), 
and is normalized as follows [8j, 



9t ^ 9L 



with the following notation: 

/ 

9l = 



= -1, 



(13) 



(14) 




(15) 



All the low-lying phonons we will consider have natural 
parity. In this case, the vertex §l possesses even T-parity. 
It is a sum of two components with spins 5 = and 
5 = 1, respectively: 

9l = gw(r)T LL0 (n,a) + g L i{r)T LL i(n,a), (16) 

where Trra stand for the usual spin-angular tensor op- 
erators [49| ■ The operators Tll and Tlli have opposite 
T-parities, hence the spin component should be the odd 
function of the excitation energy, g£' cx wj,. In this case, 
the LM amplitude in Eq. (THfl) is also the sum, 

T = T Q +T spin , (17) 

where the spin-independent LM amplitude is generated 
by the generalized EDF in Eq. ©, 



5p 2 ' 



(18) 



FIG. 2: The L-phonon creation amplitude <?i(r). 



The non-diagonal components of the matrix © , the am- 
plitudes J 7 ^^ — J 7 ^ stand for the mixed second deriva- 
tives, 



■put 



5 2 £ 
SpSv ' 



(19) 



and, finally, the amplitude is the effective pairing in- 
teraction entering the gap equation, 



(20) 



The isotopic indices in Eqs. (IT3l - |20|) are for brevity 
omitted. In the case of the surface pairing we deal, the 
"mixed" amplitude J 7 ^ is rather important for the cor- 
rect solution of the self-consistent QRPA equation (ITT 



III. PC CORRECTIONS IN MAGIC AND 
SEMI-MAGIC NUCLEI 

Let us now consider the PC corrections to Eq. ([2} in- 
duced by a L-phonon. Keeping in mind the smallness of 
the g\ parameter, we limit ourselves to the so-called g\- 
approximation. In addition, as it was mentioned in the 
Introduction, we limit ourselves with the case where the 
odd nucleon belongs to the non-superfluid sub-system. 
We will consider mainly the PC corrections for this sub- 
system, therefore the component g^ of the vector (fl"5]) 
will participate only. To shorten notation, we will omit 
for a time the upper index putting g^ — > g^ . To be more 
exact, we deal with the PC corrections to the matrix ele- 
ment Va 2 Ai = (0a 2 1 V (Ml)(j)\ 1 ) which is the second order 
variation in the phonon field g^ , see Fig. 2, 

S (2) V x . lXl = (6®cPx 2 ,Vc{> Xl ) + (cp X2 ,VSM<!> Xl ) 

with the obvious notation. In the cases where it can 
not lead to misleading, we shorten the notation below as 
follows — > 6. 



A. PC corrections included into the model 

We separate those terms of Eq. (|2"Tj) which may behave 
in non-regular way changing significantly depending on 
the nucleus and the state under considerations. As it 
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A, 




FIG. 3: A correction to the "end" 
L-phonon creation amplitude Ql(jc). 
the phonon D-function. 



The open blob is the 
The wavy line denotes 




+ 



FIG. 4: PC corrections to the mass operator, 
blob denotes the "tadpole" term. 




The dashed 



will be argued in the next subsection, this does concern 
the three first terms of Eq. (|21"j) . The last two terms 
disposed in the third line are regular and can be included 
into the TFFS parameters. As to the forth term, S^V, it 
contains both, the regular and non-regular contributions. 

The terms in the first line of (f2"Tj) . the second order 
variations of the single-particle wave functions may be 
named as the "end correction" , keeping in mind the ends 
of the diagram of Fig. 1 for the effective field. The main, 
pole part of the end correction corresponding to the sec- 
ond of these two terms is illustrated in Fig. 3. The cor- 
responding tadpole term is shown in Fig. 4. The phonon 
D-function appears here and below after connecting two 
wavy ends of Fig. 2 which corresponds to averaging of 
the product of two boson (phonon) operators B + B over 
the ground state of the nucleus without phonons. 

For the matrix element V\ 2 \ 1 of the effective field, the 
sum of these corrections corresponds to the following for- 
mula: 



sv; 



end 
A2 Ai 



— ^ V\ 2 \G\(e\ 1 )(5SaA 1 (eAj ) 
A 

^ <5£a 2 a(£A 2 )Ga(£a 2 )^AAi- 



(22) 



Here ST, means the PC correction to the mass operator 
T displayed in Fig. 4. In addition to the pole diagram 
of Fig. 4, we showed also for completeness the non-pole 
"tadpole" one Q- As we will see, this term ST tad does 



practically not contribute to the problem under consid- 
eration. 

The pole diagram in Fig. 4 corresponds to the follow- 
ing expression, 

/d 
■^~Y1 < M9lm\ x >< %lm|Ai > 



AM 



x D L (u)G x (e-u>), 



(23) 



where Dl(uj) is the phonon D-function. After an inte- 
gration, one obtains 



(e) = < X 2 \gj; M \X >< \\g L M\h > 



A M 



Tlx 



l-n\ 



£ + UJ L - e X E-UJ L - £x 



(24) 



where ujl is the excitation energy of the L-phonon and 
nx = (0, 1) stands for the occupation numbers. 
After substitution of (l24l) into (1221) one obtains: 



Va 2 a < M9lm\ x> >< x '\9lm\M > 



A' AM 



EXx + ex- 



E 

XX' M 



«A' 



£Ai — £x 
1 - n x > 

£Ai - - £a 

1 - nx 



£x 2 + — £X' £a 2 — Wi — Ex> 
< x 2\91m\ x ' >< A'|5z,m|A > Vxx t 



£a 2 



£a 



(25) 



In this equation, the term with A = Ai in the first 
sum and the one with A = A2 in the second sum are 
singular. This singularity is removed with the standard 
renormalization [7| of the single particle wave functions 
<P\i -> -JZ^ipx^ , <y5A 2 -> y/Z^ 2 ~ip\ 2 , where 



Zx = 1 - 



dVxx(e) 



ds 



(26) 



is the residue of the Green function at the pole e — ex- 

To avoid misleading, we note that below we will for 
brevity name "Z-factor" not the total quantity but the 
PC contribution only. Correspondingly, we insert into 
Eq. (H| the PC term ST only. In the TFFS the total Z- 
factor enters into the LM amplitude J- — Z 2 T^ , the local 
charges e q = ZT^ and the mean field U(r) = ZEo(f), 
where the subscript "0" means that the energy and mo- 
mentum are taken at the Fermi surface. Here the nota- 
tion Q is used with the only change a — » Z. The average 
value of the Z-factor corresponding to the total mass op- 
erator £0 was estimated in [3] and [13, EH as Zq = 0.8. 

The rest of these sums with non-diagonal terms A 7^ 
Ai in the first case and A 7^ A2 in the second can be 
calculated directly and, as we shall see, are rather small. 
However, we retain them for completeness, and represent 
the "end correction" as the sum: 



6V< 



end 
A2 Ai 



sv, 



A2 Ai 



(^aTa*)', 



(27) 
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FIG. 5: The correction due to the induced interaction. 



where 



sv x z 2Xl = (v^^-i)tW 



(28) 



Eqs. (f27j) . (f28j) correspond to partial summation of the 
diagrams of Fig. 3, and, hence, contains higher or- 
der terms in g\. To be consistent up to the order g\, 
the Z-factors Eq. (|28[) should be expanded in terms of 
<9£ A a(£)/<9£, Z = 1 + dT, X x(e)/de, with the result 



1 / <9E A2A2 ( £ ) 



de 



<9c 



13 



A2A1 



(29) 

This contribution is important for the conservation of the 
total angular momentum. This will be shown below. 

As the term 8T, tad does not depend on the energy e it 
does not contribute to Z\ and, hence, to the main term 
representing the end effect which is given by Eq. (|28[). 
It may contribute to the second term of the sum of ([27)) 
which itself is negligibly small. 

The energy derivative of the mass operator (f24|) is as 
follows: 



dSY, 



AA ( 



de 



A' M 



1 - n\> 



{e\+UJ L - E\<) 2 {S\- U L - E\,f 



(30) 



Let us now go to the second line of Eq. (|2"TT) and begin 
with the first term containing the first variations of the 
wave functions. After connecting two wavy ends we ob- 
tain the "triangle" diagram (GGD) displayed in Fig. 5 
which describes the effect of the induced interaction .Find 
due to exchange of the L-phonon. The explicit form of 
the corresponding correction to the matrix element V\ 2 Ai 
is 



A A' M 

x G x (e Xl - uj)Gy{e\ 2 -u)D L (u). 



FIG. 6: Correction due to the direct action of the Ml-field 
to the phonon. The black blob is the magnetic moment of the 
phonon. 



^X = (-ir 2 - ro2 (_t 2 m ^)<nsv—\\i>, 

(32) 

where the short notation 1 = v\ is used and the reduced 
matrix element is equal to 

< 2 || SV GGD || I >= Y(-l) L +l*-h { 1 1} h 
x <4|| V||3><2|| g L ||4><3|| g L || 1 > 7 34 (w L ), (33) 



1 



(e 3 - £4) - (£1 - £2) 
77,4 I — 713 



"3 



Si - £3 + Wi 

1 — 714 



£2 - £4 + u L £1 - £3 - uj l £ 2 - £4 - Wi 



(34) 



where 



9l = 9l(-u) = gLo(r]Uj)T LLO (n,a) 

~g L i{r;u>)T LL i(n,a), (35) 

Let us go to the second term in the second line of Eq. 
(f2"T|) with the second order variation of the effective field 
V. This term was considered in detail in (53[. We sep- 
arate the PC correction due to the phonon contribution 
to the local charge e q in the first term of Eq. ([2]) for 
the effective field. It is represented by the diagram of 
Fig. 6 describing the contribution of the phonon mag- 
netic moment. It can be obtained from the first diagram 
of Fig. 4 by inserting the external magnetic field in the 
phonon line. In Fig. 6, the black blob means the phonon 
magnetic moment, 



ph 

Ml 



9?L, 



(36) 



= /|£ £ <*2\9l M \X>Vxx<\\g LM \\ 1 > 

(31) 



where g£ h is the phonon gyromagnetic ratio. In the 
Bohr-Mottelson model [4l|, [52|, one has g£ h BM = Z/A. 
Similar to to Eq. pip , we get 



After separating the angular variables and integrating 
over the energy, we obtain: 



[ £ J2 < A 2|.9lmI a >< %lm|Ai > 

J A M 



x G a (£a 1 - uj)D l (u})D l (u - w ). 



(37) 
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FIG. 7: The tadpole-like diagram for the contribution of the 
phonon magnetic moment. 



Again, after separating the angular variables, we ob- 
tain for the reduced matrix element, 



< 2 || SV GDD || 1 >= 



■33 



1L 



x <2|| g L || 3 > (l 3 (1) (o; i )+7 3 (2) (a; i )) , (38) 

l - m 



(ei - e 3 - uj L )(e 2 - £3 - u L ) 

m 

{ex - e 3 + u> L )(e 2 - £ 3 + wl) ' 



(39) 



£i-£2-2wl V£2-£3+^L 
1 / 1 - n 3 



1 f n% 1 — n 3 

£l-£ 3 -WL 

£l-£ 2 + 2w L V £ 2 - £3 - W L ' £i-£ 3 +CJ L 



(40) 



The second integral (|4T)|) reveals a dangerous behavior at 
(ljl — wo) — > 0. For the case of the static external held 
we deal, it occurs at ljl 0. 

Analogically to producing Fig. 6 by inserting the ex- 
ternal held to the first, pole, diagram of Fig. 4, there is 
the tadpole-like PC correction to the effective magnetic 
field which corresponds to inserting the external field to 
the phonon D-function in the second diagram on Fig. 4. 
It is displayed in Fig. 7. The dashed blob denotes the 
sum of all non-pole diagrams of phonon-particle scatter- 
ing (the phonon-particle Compton diagrams). It contains 
the integral 



I 



tad 



2vri 



D l (uj)D l (uj - wo) 



4a; j 



(41) 



Its behavior at u>l —> is just the same as of the integral 
J 3 , Eq. (|4"0")l . This makes it reasonable to suppose that 
their sum is regular at (w/, — loq) — >• 0. 






FIG. 8: 
nuclei. 



Diagrams for the phonon magnetic moment in magic 



The phonon gyromagnetic ratio <?£ h for magic nuclei is 
obtained from the diagrams shown in Fig. 8. For non- 
magic nuclei, the corresponding set of diagrams is much 
more complicated, see Ref. [H, HH for the case of the 
quadrupole moment, and in this case we limit ourselves 
to the BM model. 

Thus, the sum of the PC corrections to the effective 
M 1-field is as follows: 



SV 



SV Z + SVggd + SV, 



(i) 



GDD 



SV' 



end 



SV 



(2) 



GDD 



SV 



non— pole 



(42) 



where SVq^' d 2 ^ are the terms of Eq. (|3"T|) with the inte- 
grals /W'( 2 ) and SV non - po \ c denotes the sum of all non- 
pole PC corrections to the magnetic effective field. Up 
to now, the method of a consistent calculation of this 
quantity is implemented for magic nuclei mainly 0, l42j . 
For non-magic nuclei, the formalism was developed in 
[HI but the only one realization was carried out [33[ . As 
we shall see below, the magnitude of the first two terms 
of the sum (|42|) is, as a rule, significantly larger than 

that of SVqq^ terms. However, SV Z and SVggd terms 

are always of the opposite sign and strongly compensate 

each other. Therefore, the term SVq^ d is sometimes not 

negligible and should be taken into account. In the case 

when the external field is Vq — j, where j = 1 + l/2er 

is the total angular momentum, due to the conservation 

of j condition, the total PC correction (|42l) should be 

zero. As we will discuss in the next section, the sum of 

the first four terms of (|4"2"]) is equal to zero. This means 
(2) 

that the terms SV G ^ D and <5V^ lon _ po ic should compen- 
sate each other. We will suppose that a similar relation 
is approximately valid for any Ml-field and will use the 
model where the terms in the first line of Eq. (|4"2"j) are 
taken into account only and those on the second line are 
neglected. In fact, we carry out the summation of the 
end corrections, similar to that in Eq. ()28|) , and the final 
ansatz we use is as follows: 



V\ 2 \ t — \J Z\ 2 Z\ X 

x (V + SV GGD + SVV D +Sv: nd ) .(43) 

\ / A2A1 

Just as Eq. (J2HJ) vs Eqs. fl23]) and (US]), the ansatz 
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differs from the prescription of Eq. (|42|) in terms higher 



ln 9l- 



In conclusion of this section, we write down the ex- 
pression for the L-phonon magnetic moment /i^ in magic 
nuclei corresponding to Fig. 8: 



/i L =^(-l) 

123 



L + l 





L L\ 




L 




c 


L —L J 


\33 


.12 





jGGG 
xi 123 



x < 1 || V(M1) || 2 > [< 1 || g L || 3 >< 3 
K) 



9L || 2 > 



where 



<1|| g L || 3><3 1| g L ||2> /gf G (-u L )] ,(44) 
- n 2 )(l - n 3 )-(l-ni)n 2 n3 



£ 2 - ei 



£l - £3 - WL 



ni(l - n 2 )rt3 - (1 - ni)n 2 (l - 713) 



£2 - £3 - u L 
nin 2 (l - n 3 ) - (1 - ni)(l - n 2 )»3 
(£1 - £3 - Wi)(e 2 - £3 - wl) 



(45) 



For the general case of luq 7^ 0, uj' l = ujl + ujq =/= ljl 
formulae analogous to Eqs. (|44|) and (145)) for the E2 
symmetry have been obtained in Ref. [551 ] . 



B. PC corrections we omit 

The diagram in Fig. 6 corresponding the variation of 
the local charge e q in the effective held V is the only 
one example of g\ PC corrections to the effective held 
S( 2 >V. The main part of such corrections appears from 
variations of the integral term of Eq. (|2l). These terms 
were examined in detail in Refs. [43|, Hj, [56| . The term 
corresponding to the variation of the particle-hole propa- 
gator, 6^A(1) = S^GS^G is displayed in Fig. 9. The 
wavy line for the phonon D-function, just as above, ap- 
pears after connecting the phonon ends of Fig. 2. This 
diagram, together with that in Fig. 5, can be interpreted 
as the substitution of the sum J rs P in — y jrspm _|_ jr ;nd as 
the effective interaction in Eq. ^ with different order 
of these operators. Another variation of the particle-hole 
propagator, S^A(2) = 2S^GG, is displayed in Fig. 10. 
There are also more complicated diagrams for S^V con- 
taining mixed variations ^^S 1 - 1 "' A, see (53|. 

Let us discuss the reasons why we exclude these two 
diagrams from our model although they are of the same 
g\ order. In fact, as it is seen from all the equations 
of the above subsection, g\ appears always with some 
energy dependent denominators, and the characteristic 
dimensionlcss combination is the ratio of 



(e x - £a' ± w l ) 2 



(46) 



Typical value of the matrix element (<?z)aA' — 1 MeV 
whereas the energy differences in the denominator of (|46[) 
are of the order of e F /A 1/3 ~ 7 -j- 10 MeV. It should 
be reminded that we deal with "magic" sub-system of a 




FIG. 9: The diagram alternative to the one in Fig. 5. 




FIG. 10: A PC correction to the particle-hole propagator in 
Eq. ©. 



semi- magic nucleus where pairing is absent. Therefore, 
as a rule, (3l is really a small parameter and the il g\ 
expansion" is justified, each separate term is small and 
only a large sum of such contributions can be important. 
Such sums do practically not depend on the nucleus un- 
der consideration and on the single- particle state of the 
odd nucleon. As simple estimations show, this is true 
for diagrams in Fig. 9 and Fig. 10. According to the 
strategy of our model, we must consider them as contri- 
butions to the LM amplitude .F spin and the local charge 
e q and not calculate explicitly. On the other hand, as 
we shall see below, the PC corrections considered in the 
previous Subsection may depend, sometimes strongly, on 
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i M1 




FIG. 11: An example of more complicated "mixed" diagrams 
corresponding to the third line of Eq. <|21[) . 



the state Ao of the odd nucleon with rather big contribu- 
tion of a separate term with small energy denominator 
in Pl. Such PC corrections depend significantly on the 
nucleus under consideration and on the state Ao of the 
odd nucleon. They should be calculated explicitly. 

Fig. 11 represents an example of the "mixed" PC cor- 
rections from the third line of Eq. (12T1) . All arguments to 
exclude the two diagrams from our model remain valid 
for them too. 



IV. PHONON MAGNETIC MOMENTS IN 208 PB 

In this double-magic nucleus several low-lying phonons 
are known with varying degree of collectivity. The low- 
est one 3~ possesses the highest collectivity whereas the 
strength of the 5~ mode is shared between two states. 
Before starting the study of PC effects in 208 Pb we ana- 
lyze the accuracy of describing the phonons themselves. 
Their characteristics are presented in Table I. We see 
that the most collective 3~ level is described sufficiently 
well, 5j~ and ones, a little worse. As to the 2^~ state 
and other states of positive parity, they are not strongly 
collective as there is no low-energy particle-hole con- 
figurations of positive parity except spin-orbit doublets 
(h^n hg/2) for protons and (« 1 ~ 3 1 /2*ii/2) f° r neutrons. In 
such a situation, ujl values in the QRPA solution are 
shifted only little from the nearest particle-hole excita- 
tion energy, e p (h g /2)-ep(h 11 / 2 ) or e n (i n / 2 )-Sn(ii3/2) in 
our case. Single-particle levels in 208 Pb generated with 
the DF3-a functional we use are displayed in Fig. 12 
for neutrons and Fig. 13 for protons. Comparison is 
made with the experimental spectra and those obtained 
by us with the one of the best Skyrme functionals HFB- 
17 0. We see that our spectra reasonably agree with the 
data, on average, better than the HFB-17 one. However, 
some inaccuracy takes place for levels under considera- 
tion which is partially responsible for the shift of the 2^ 
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FIG. 12: Neutron single-particle levels in 208 Pb. 



TABLE I: Characteristics of the low-lying phonons in 

2r 2L, 



u) L (MeV) and B(EL, up)(e mi 



B(EL) th B(EL) C 



3r 

■A 

4 



2.684 
3.353 
3.787 
4.747 
5.004 
4.716 



7.093 x 10 J 



2.615 
3.198 
3.708 
4.086 
4.928 
4.324 

5.367 4.911(?) 8.462 x 10' 
4.735 - 6.082 x 10 
5.429 - 1.744 x 10 1 



3.003 x 10 8 
1.785 x 10 8 
1.886 x 10 s 
1.148 x 10 3 
3.007 x 10 6 



6.12 x 10 5 
4.47 x 10 8 
2.41 x 10 8 
3.18 x 10 3 



level up to 660 keV. Another reason for this discrepancy 
is evidently not taking into account the spin-orbit LM 
amplitude in Eq. (|T3)) for the gi, vertex. We shall see 
below that the PC corrections in magic nuclei originate 
mainly from the 3~ phonon, therefore some inaccuracy in 
describing higher collective states does not lead to serious 
errors. 

Let us now calculate magnetic moments and corre- 
sponding gyromagnetic ratios g£ h of the i-phonons in 
208 Pb. It is given by the diagram of Fig. 8, Eq. (pHl) . 
We need the g£ h values to evaluate the contribution of 



the phonon magnetic moment term, Fig. 6 and Eq. 
In addition, this calculation elucidates to a large extent 
the structure of the phonons under consideration. The 
results are given in Table II. We showed separately the 
j- and s-components, according to Eq. ([TJ , with obvious 
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TABLE II: Magnetic moments (in /ijv units) of phonons in 208 Pb. 



LT 




W 






M 

r i 1 


Lin 
f-'P 




M 

Ut 

' L 


LIT 




3" 


-0.074 


-0.039 


-0.113 


1.566 


0.058 


1.492 


1.492 


0.019 


1.511 


0.463 


sr 


-0.027 


-0.018 


-0.046 


4.733 


0.278 


5.011 


4.705 


0.260 


4.965 


0.993 


5,7 
2 


-0.215 


-0.478 


-0.693 


0.853 


-0.123 


0.730 


0.638 


-0.600 


0.037 


0.008 


'2t 


-0.027 


0.000 


-0.027 


1.536 


0.493 


2.029 


1.509 


0.492 


2.002 


1.001 


2+ 


-0.027 


0.004 


-0.022 


1.541 


0.406 


1.947 


1.514 


0.411 


1.925 


0.962 




-0.009 


-0.010 


-0.018 


4.017 


0.449 


4.466 


4.008 


0.440 


4.448 


1.112 




-0.112 


-0.232 


-0.343 


1.822 


0.276 


2.098 


1.711 


0.044 


1.755 


0.439 




-0.005 


-0.004 


-0.009 


6.172 


0.294 


6.466 


6.167 


0.290 


6.457 


1.076 


6 2 + 


-0.075 


-0.147 


-0.222 


4.765 


0.092 


4.857 


4.690 


-0.054 


4.636 


0.773 
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FIG. 13: Proton single-particle levels in 208 Pb. 



substitution of 1 = j — s. The subscripts n, p refer to the 
neutron and proton subsystems whereas L corresponds 
to their sum, e.g. yq) = ^ + [ip\ To understand 
better the nature of phonons we investigate it is worth 
comparing the g£ h values obtained with the BM model 
prediction g£ h BM = Z/A = 0.394. We see that only for 
the 3~ and 4^"-states our values are rather close to the 
BM one whereas in other cases there is nothing in com- 
mon with these two theory predictions. Note that in 

( s) 

the BM model the spin-component /j} L ' is absent. If we 
neglect in Eq. (pH)) the spin term of the effective field 
V(M1) and put = 0, i.e. take V(Ml) = j, we obtain 
the BM value of g£ h for all the states under considera- 



g(3"), MeV 




FIG. 14: Components of the vertex gL , V = 3 _ , in 208 Pb. 



tion. Thus, the microscopic value of the gyromagnetic 
ratio deviates from the classic model prediction due to 
the spin term and non-zero value of Q . 

To testify the formulas above and estimate the accu- 
racy of the calculations, it is instructive to apply Eqs. 
(gU, (gSJ) to the spurious phonon 17 = 1~, uj 1 - = 0. In 
this case, the term gn in the sum of Eq. (j 161) vanishes, 
whereas the term gio(u>) is singular at ui — [7j, 



010 M = 



1 



dU 



(47) 



where U(r) is the central part of the mean- field poten- 
tial generated by the energy functional (j9|) and B\ — 
3m/(47rA) is the BM mass coefficient 01 . Eq. (g7J) fol- 
lows from the exact TFFS self-consistency relation Q 
with some simplifications and neglecting the spin-orbit 
terms. The singularity in all the above expressions con- 
taining g\ is compensated as the corresponding integrals 
of the Green functions are proportional to w. This ap- 
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creation amplitudes gL are displayed in Figs. 11-13 for 
the states 3~, 5^ and 2^, respectively. We see that all of 
them contain the BM-like (oc dU/dr) surface maxima of 
the spin-zero components which are significantly larger 
than components with S = 1. However, for 5^ and 
states they are not negligible. Moreover, they possess 
maxima at r ~ 6 fm where the wave functions of the 
single-particle states in vicinity of Fermi level have their 
maxima too. This is true for couples of the wave func- 
tions (<fi\4>\i) which correspond to the term with a small 
energy denominator in Eq. (|45p which always exists for 
the phonon states under discussion. At the same time 
the maxima of the main, S = 0, components are shifted 
to the right at ~ 1 -j- 2 fm due to the coherent contribu- 
tion of a lot of single-particle states distant from Fermi 
surface. However, in this region the wave functions under 
discussion began vanish. As the result, for these states 
the contribution of the S = 1 to the phonon magnetic 
moment is often rather big. 



FIG. 15: Components of the vertex g L , L n =5^, in 208 Pb. 



g(5"), MeV 
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r, fm 

FIG. 16: Components of the vertex g L , V = 2+, in 208 Pb. 



proximation for gi violates a little the normalization re- 
lation (fl"4")l leading to the value -1.074 instead of -1. If 
we correct normalization, we obtain for the magnetic mo- 
ment of the spurious l~-phonon we obtain /i(l~) = 0.403 
in very good agreement with the BM gyromagnetic ratio 
0.394. This calculation confirms the self-consistency of 
the scheme developed above. 

In our calculations, the spin component is negligible for 
the 3 _ -state only. It confirms that this state, indeed, is 
most similar to the BM surface vibrations. The phonon 



V. J-CONSERVATION CONDITION 

As the total angular momentum j = 1 + l/2er is the 
integral of motion, the identity 

V(j)=j (48) 

should be fulfilled. Within the standard TFFS, it leads 
to relations for the local charges e" ,p ||. Obviously, when 
the PC corrections are included, the sum of all of them 
discussed above should vanish in the case of Vq — j. As it 
was mentioned in Sect. 3, the sum of the first four terms 
of (|4"2"]) is equal to zero for any value of L. This can be 
shown analytically and by numerical calculations as well. 
This means that the identity 

<SKon- P oic(j) = ~5V£ ] DD {S) (49) 

should be fulfilled to guarantee the identity (|4"5|) . We 
suppose that the similar relation, 

^Kon-poic(Ml) = -5V£ ] DD {M1), (50) 

is valid for the complete Ml-field. This relation is the 
main ingredient of the model we use. 

It is worth noting that the ansatz (|43[) violates the 
j-conservation law as it takes into account higher order 
terms in g\ in an non-consistent way. However, this vio- 
lation turned out to be rather small. 



VI. MAGNETIC MOMENTS OF THE ODD 
NEIGHBORS OF 208 PB 

As it was mentioned above the first two terms of 
Eq. (|4"2")l dominate. Let us begin from the "end cor- 
rection" which is simpler for the analysis. In Table III 
we show contributions of three phonons to the Z-factors 
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TABLE III: Contributions of different low-lying phonons to 
the Z-factors of single-particle states in 208 Pb. 



A 


z(3r) z(5r) ^(2+) 




z' 


Ztot 


3^5/2 





.963 





.989 


0. 


.967 


0. 


.944 





.869 


i -n 
1J15/2 


0. 


.657 


0. 


.985 


0. 


.989 


0. 


.974 


0. 


.623 


1*11/2 


U 


07/1 


U. 


yyo 


U. 


yyu 


0. 


.981 


U. 




2<?9/2 


n 


Q1 8 


n 
u. 




n 
u. 




0. 


.975 


n 
u. 


881 
001 


3pr; 2 





.963 


0. 


.995 


0. 


.986 


0. 


.979 





.925 


3^3/2 





.960 


0. 


.992 


0. 


.984 


0. 


.971 





.910 


2/5-/2 





.964 


0. 


.997 


0. 


.985 


0. 


.974 


0. 


.922 


1 A n 
- Lt 13/2 





.943 


0. 


.990 


0. 


.990 





.979 





.904 


2/ 5 P / 2 





.953 


0. 


.995 


0. 


.941 


0. 


.907 


0. 


.810 


li p , 

13/2 





.747 





.992 





.993 





.975 





.718 


2/r P / 2 





.898 





.996 





.991 





.971 





.861 







.978 


0. 


.999 


0. 


.996 





.987 





.960 


O -P 

°*l/2 





.959 


0. 


.996 


0. 


.992 


0. 


.979 


0. 


.928 


ZU 3/2 





.964 


0. 


.998 


0. 


.992 





.982 


0. 


.937 


1 ' 1 ll/2 





.951 





.996 





.995 


0. 


.982 





.926 


2^2 





.758 





.996 





.986 





.968 





.721 



of different single-particle states in the odd neighbors of 
208 Pb, i.e. in 2 07,209p b i sot0 pes for neutrons and 207 T1, 
209 Bi nuclei for protons. Here, we used the notation 
Ztot — Yli—i Zi and Z' = Z3 Yli=5 Z%, where the index i 
numerates all the collective states listed in Table I. Thus, 
the states 37/, 57/ and 2^~ are excluded in Z' . 

As it was stated in the Introduction, the contribution 
of the 3~ state to the PC corrections dominates in magic 
nuclei due to its strong collectivity. Table III confirms 
this statement on average, although often the the sum 
of the contributions of all other phonons to the differ- 
ence of (1 — Z\) is comparable with the one of the 3~ 
phonon alone. From the other phonons, the 2^~ phonon 
contributes most to the Z-factor although it is less col- 
lective than 5j~, 5.7 phonons. The reason becomes clear if 
one looks to Eq. (|30l) which determines the phonon con- 
tribution to the difference of (1 — Z\), see also Eq. (|46]l 
which gives an estimate of a separate contribution to dif- 
ferent PC correction. It is obvious that the situation is 
on average more preferable for phonons of positive par- 
ity. Indeed, in this case the states A, A' in ([35]) are of the 
same parity and it is easier to find those single-particle 
states with a "small denominator". That is why there 
are two cases, the 3dgy 2 and 2/^ 2 states, where the con- 
tribution of the 2^ phonon to the difference of (1 — Z\) is 
comparable with the one of the 3~ phonon, although the 
matrix elements of 172 are typically three times smaller 
than those of (73. Typical value of the Z- factor in 208 Pb 
is Ztot = 0.8 -i- 0.9. In this case, the perturbation theory 
in g 2 should work well. However, there are one neutron 
state, lj™5/ 2 : an d two proton ones, ljf 3 y 2 an d for 
which we have Z to t = 0.7 -r 0.6 and application of the 
perturbation theory is questionable. 

Table IV contains various PC corrections to magnetic 



moments of odd neighbors of 208 Pb. In agreement with 
the above ansatz, the total value of Sfi given in the last 
column is the sum of the terms in the first line of Eq. (|42p . 
The term Sfi GDD is presented just for information. In the 
upper part of the table five cases are presented for which 
experimental data exist. We see that for all of them the 
PC correction to fi value is negligible. The reason is the 
strong cancelation of the first two corrections, Sfi z and 

SfiQGD- The term 6ivq DD is usually significantly smaller 
than these two main ones. The term 5/j,' end is defined 
in accordance with Eq. (|2"T)) . This correction is always 
negligible and is included just for completeness. 

In the two bottom lines of Table IV we presented two 
examples of strong PC corrections for two states which 
magnetic moments are not known. In both cases, the 
PC correction to the Z-factor is rather big, see Table III, 
because of the existence of a term with small energy de- 
nominator, see Eq. (|46[) . for the 3~-phonon. Due to the 

same reason, the term 5^q^ dd is also rather big. In the 
case of the lj™ 5 / 2 state, Sfi and S^ggd terms again can- 
cel each other, but the 5^q DD term is comparable with 
the first two terms and leads to a big total correction. In 
the 1^3/3 case, the cancelation is not so strong as usual, 

but the 5^q DD term is again important. 

Table V contains the final results for magnetic mo- 
ments, 

Mth = Mo + Sfi, (51) 

where the first term is the TFFS prediction for the mag- 
netic moment, i.e. the solution of Eq. and Sfi is 
the sum of all the PC corrections considered above. We 
see that they can not help in solving the long-standing 
problem of the theoretical description of the ground state 
magnetic moment of the 209 Bi nucleus. Evidently, in this 
case the usual TFFS form © for the spin-dependent LM 
amplitudes is not sufficient and the spin-orbit terms and 
tensor ones in cross channel should be included. This 
is rather difficult in the coordinate space method we 
use for solving Eq. @. The direct solution in the A- 
representation in a large basis would be more adequate 
in this case. Importance of the spin-orbit force for this 
problem was found in Ref. (57j , 

In conclusion of this Section, we repeat that our model 
for PC contributions to magnetic moments of odd neigh- 
bors of 208 Pb leads to very small corrections. A similar 
result was obtained previously by Tselyaev [43| who con- 
sidered mainly diagrams we omit, displayed in Fig. 9 and 
Fig. 10. The reason was in strong cancelation of these 
diagrams similar to cancelation of the diagrams in Fig. 3 
and Fig. 5 within our model. 

VII. MAGNETIC MOMENTS OF SEMI-MAGIC 
NUCLEI 

Let us go to semi-magic nuclei. As it was mentioned 
in the Introduction, we deal with such odd nuclei where 
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TABLE IV: PC contributions to the magnetic moments (in /j,n) of odd neighbors of 208 Pb. The notation is explained in the 
text. 



nucleus 


A 






5n z 


Suggd 


r^GDD 


( s Pahn) 




5/ipc 


207p b 


3PV3 


0.474 


3r 

all 


-0.018 
-0.037 


0.015 
0.030 


-0.2E-3 
0.006 


0.018 
0.031 


0.5E-3 

0.002 


-0.003 
0.001 


207p b 


2/ 5 72 


0.720 


sr 

all 


-0.027 
-0.058 


0.019 
0.046 


-0.001 
0.050 


0.071 
0.140 


-0.002 
-0.002 


-0.011 
0.036 


209p b 


2^9/ 2 


-1.337 


sr 

all 


0.120 
0.174 


-0.111 
-0.154 


-0.051 
-0.013 


0.166 
0.157 


0.007 
0.008 


-0.035 
0.015 


207r|-ij 


Q e p 
° b l/2 


1.857 


3i 
all 


-0.082 
-0.137 


0.089 
0.130 


-0.004 
-0.005 


0.005 
0.001 


-0.004 
0.035 


-0.001 
0.023 


209 R - 

1)1 


1 h p 

1,i 9/2 


o.uyi 


all 


-u.uoo 
-0.146 


U.UJl 

0.096 


n m a 
0.033 


0.041 


u.uuz 
0.002 


-u.uio 
-0.015 


209p b 


lJlS/2 


-1.315 


Si- 
all 


0.688 
0.754 


-0.679 
-0.738 


0.696 
0.780 


-0.237 
-0.249 


-0.001 
-0.001 


0.703 
0.794 


209 Bi 


li p 

"l3/2 


8.071 


sr 

all 


-2.730 
-3.020 


1.530 
1.730 


0.413 
0.460 


-0.162 
-0.178 


0.3E-3 
0.001 


-0.785 
-0.824 



TABLE V: Magnetic moments (in /ijv units) of the odd neigh- 
bors of 208 Pb. 



nucleus 



/'() 



dfj, 



A^cxp 



207p b 
207j 
209j 
209j 

207pj 

209 Bi 
209 Bi 



7 Pb 

3 pb 
'Pb 



3?! 

2/, 



/2 
— n 
5/2 

2<?9/2 

ljlS/2 
o -P 
°*l/2 

lh p 
li p 

±l 13/2 



0.474 0.001 

0.720 0.036 

-1.337 0.015 

-1.315 0.794 

1.857 0.023 
3.691 



0.592585(9) 

0.80(3) 
-1.4735(16) 

1.876(5) 

-0.015 3.676 4.1106(2) 



0.475 
0.756 
-1.322 
-0.521 
1.880 



8.071 



-0.824 



7.247 



the odd nucleon belongs to the non-superfluid subsystem. 
Otherwise, all equations of Sect. 3 should be generalized 
by including the superfluidity. In particular, we consider 
the odd-proton neighbors of the even lead nuclei, i.e. Tl 
and Bi isotopes, and the same for the tin-core region, 
where we analyze In and Sb isotopes. However, the neu- 
tron sub-system which is now superfluid and we have 
to apply the vector notation of Eq. (|15|) using the ini- 
tial notation gy^ for the normal component of <?l and 

g L ' for the anomalous ones. All the equations of Sec- 
tion 3 remain valid with the only exception of Eqs. (l4"4"l) . 

corresponding to Fig. 8, for the phonon magnetic 
moment or the gyromagnetic ration g£ h . Indeed, now 
the neutron triangle should be included also, with taking 
into account the superfluidity effects. Therefore we use 
the BM model prediction for g£ h in this Section. This is 
justified as far as we consider only such phonons which 
are very collective. 

In all even-even spherical non-magic nuclei there exists 



TABLE VI: Characteristics of the low-lying 2 ] f -phonons in 
even Pb isotopes, uj 2 (MeV) and B(_E2,up) x 10 4 (e 2 fm 4 ). 



A 



U>2 



B(E2) th B(E2) C 



188 1.028 
190 0.930 
192 0.849 
194 0.792 
196 0.764 
198 0.762 
200 0.789 
202 0.823 
204 0.882 
206 0.945 
208 4.747 
210 1.346 
212 1.443 
214 1.125 



0.723 
0.733 
0.853 
0.965 
1.049 
1.063 
1.026 
0.960 
0.899 
0.803 
4.086 
0.799 
0.804 
0.836 



0.551 
0.617 
0.634 
0.646 
0.627 
0.624 
0.479 
0.373 
0.250 
0.130 
0.189 
0.036 
0.131 
0.161 



0.162 (0.004) 
0.100 (0.002) 
0.30 (0.03) 
0.051 (0.015) 



a very collective low-lying 2 J"-state with an excitation en- 
ergy u>2+ — 1 MeV. In fact, it is the only state which 
should be taken into account for the PC corrections to 
magnetic moments. For the odd-proton neighbors of nu- 
clei soo-^ph we included also the 3^-state, the next in 
the range of collectivity, and found that it gives only sev- 
eral % of the main PC correction of the 2f -state. There- 
fore in other cases we limit ourselves with this PC cor- 
rection only. In Table IV we present excitation energies 
and B(E2) values of the even Pb isotopes. 

These results were reported previously in [l5j . Exper- 
imental data on B(E2) values are not complete whereas 
the w 2 + value is known for all isotopes in the chain we 
consider. As the analysis in the previous Section showed 
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TABLE VII: Z-factor values of two proton levels close to the 
Fermi surface due to the 2f and 3j~-phonons in the lead chain. 

A Z(3s 1/2 )[2+] Z(3s 1/2 )[3-} Z(lfy a )[2+] Z(l/i 9/2 )[3-] 



188 


0.830 




0.659 




190 


0.822 




0.614 




192 


0.834 




0.605 




194 


0.719 




0.542 




196 


0.824 




0.547 




198 


0.810 




0.526 




200 


0.755 


0.966 


0.619 


0.980 


202 


0.743 


0.967 


0.628 


0.981 


204 


0.856 


0.968 


0.786 


0.982 


206 


0.915 


0.965 


0.880 


0.981 



g <0) (2 + ), MeV 
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- P, s=o 
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FIG. 17: Components of the normal amplitude g 2 in Pb. 



the latter is of primary importance for calculating the PC 
corrections. We see that the excitation energies on aver- 
age are in reasonable agreement with the data. However, 
the theory does not reproduce the experimental trend. 
Specifically, if we exclude the magic 208 Pb nucleus, the 
experimental uj 2 + values have a maximum, ~ 1 MeV, 
in the middle of the chain, A ~ 200, and decrease to 
~ 0.7 -T- 0.8 MeV at both ends. The theoretical predic- 
tions have opposite behavior. They show in the middle 
a minimum, ~ 0.8 MeV, and increase up to ~ 1 MeV in 
the left end and up to ~ 1.4 MeV in the right one. There- 
fore we may expect that our predictions for S/ipc values 
will be too high in the middle of the chain and too low 
in the both ends. However, as the estimation (|4"oT) shows, 
the position of single-particle levels in the Fermi surface 
vicinity are often as important as the u>l value. In the 
very right end of the chain the difference of uif 1 — 
reaches so large value as 0.6 MeV. In view of such de- 



g (1 ' 2) (2 + ), MeV 




J I I I I I I I I I I I I I I I I I I I L 



2 4 6 8 10 

r, fm 

FIG. 18: Components of the neutron anomalous amplitudes 
and g { 2 2) in 200 Pb. 



ficiency, we may wait most errors of our predictions for 
PC corrections in this part of the chain. 

In Fig. 17 we display the components of the normal 
amplitude gf^ in 200 Pb which is in the middle of the 
lead chain. For graphical convenience the spin-one com- 
ponents are multiplied by the factor of 5. Remind that, in 
fact, we deal with the proton vertices only. We see that, 
just as for the 3~-state in 208 Pb, the spin-one component 
is negligible and can be neglected. For completeness we 
present also in Fig. 18 the neutron anomalous ampli- 
tudes g^\ spin-zero components only. We see that 
the first of them is rather big, comparable to the normal 

vertex gf\ Note that all the components of the vector 
gL are normalized with the same normalization condition 
(JH and can be compared directly. 

To estimate the scale of the PC corrections, it is in- 
structive to analyze the Z-factors of the states we con- 
sider. We deal with the odd Tl isotopes with 3si/ 2 state 
for the odd proton and the odd Bi isotopes with I39/2 
state for the odd proton. Corresponding values of the 
Z-factors induced by the 2f -state are given in Table VII. 
For lead cores A — 200 4- 206 we showed also the Z- factor 
induced by the 3]~ -state. We see that it is rather close 
to unit, so that Z tot = Z(2 + )Z(3") ~ Z(2+). This is 
a signal that all the PC corrections due to the 3^~-state 
are now very small and those from the 2^~-state could be 
taken into account only. 

We see that for the major part of the lead isotopes the 
Z-factors deviate from unit significantly, Z ~ 0.5 -j- 06. 
This occurs because in the sum of Eq. (f3"0"|) which deter- 
mines the Z-factor there is one or two "dangerous" terms 
with small energy denominators for which the perturba- 
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tion parameter (|46[) is not small. In such a situation, 
the perturbational in g\ theory becomes doubtful and 
one has to consider higher order corrections. There is a 
simple way to take into account some of them. Let us 
separate the PC correction to the effective field, the first 
line of Eq. (|42|) . into two components, 



SV = 5V Z + SV, (52) 

and the PC correction to the magnetic moment, 

Sfi = Sfi z + Sfi'. (53) 

Instead of the direct perturbation formula for the mag- 
netic moment (fBTj) we will use also the following ansatz: 



Ath = Zx(fi + 6fi'). 



(54) 



The physical meaning of this recipe is very simple. It 
corresponds to introducing the "end correction" similar 
to that of Fig. 3 not only to the effective field V(M1) 
itself, but also to triangle diagrams of Fig. 5 and Fig. 6 
describing other types of the PC corrections. Of course, 
this is just an intuitive ansatz and a consistent way to 
consider the dangerous terms is necessary. Up to now 
it is not developed and we present in Fig. 19 both pre- 
dictions for Tl isotopes of Eq. (jSTI) corresponding to the 
g^-approximation and of the ansatz (I54[) where also some 
higher in g\ terms are considered. In this case, the abso- 
lute value of the PC correction of Eq. (1541 is reduced, al- 
though the difference between these two results is rather 
small. 

We see that for all isotopes of this chain except 207 T1, 
i.e. all non-magic ones, the sign of the PC correction is 
negative which always helps to reproduce the data. The 
experimental fi value undergoes a sharp jump from 207 T1 
to 205 T1 remaining almost constant for lighter isotopes. 
The PC correction describes such a behavior only quali- 
tatively. Its absolute value grows smoothly till 199 > 201 T1 
and becomes again rather small for lighter isotopes. Ta- 
ble VIII shows separate components of Sfi, with the same 
notation as in Table V. Again, as it took place for the 
odd neighbors of the magic 208 Pb, the resulting Sfi value 
is much smaller than the terms Sfi z and Suggd sepa- 
rately which have opposite signs and compensate each 
other strongly. Other components of Sfi again are less 
important. However, the sign of the 5fi^) DD term is the 
same as of S/iggd and "helps" to diminish the total Sfi 
value. 

Similar results for Bi isotopes are displayed in Fig. 
20 and separate terms of Sfi are presented in Table IX. 
Now the cancelation of different PC corrections is even 
stronger than for the Tl isotopes, and the \5fi\ value is 
usually 10 or 15 times less than |<5/u z |. Contrary to the Tl 
case, now the sign of the PC correction looks incorrect. 
However, the main reason for the contradiction between 
the theory and experiment is evidently incorrect value of 
fiQ. Indeed, in 209 Bi where the PC correction is negli- 
gible, see Section 5, we have fi c ^ p — A*o — 0.5 fiN- As it 



M. Mi 




187 191 195 199 203 207 211 

A 

FIG. 19: Magnetic moments of Tl isotopes. 
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FIG. 20: (Magnetic moments of Bi isotopes. 



was discussed above, the main reason for this deviation is 
supposedly absence of the spin-orbit and additional ten- 
sor terms in the spin-dependent LM interaction ([5]). In 
any case, the final conclusion about the sign of the PC 
corrections could be made only after solving the problem 
of 209 Bi. 

Let us go to the tin region. Again we consider the 
proton-odd neighbors of even Sn isotopes, the odd In 
and Sb isotopes, with non-superfluid proton sub-systems. 
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TABLE VIII: PC contributions to the magnetic moments (in (j,n) of odd Tl isotopes, Ao = 3s 1 ? 2 . The notation is explained in 
the text nearby Table IV. 



A a 








17 


U fJj 


O^GGD 


f^GDD 


(Su {2) ) 
("h-gdd) 


tin' 




Mth 


Ath 


207 +1.876(5) 


1 


857 




all 


-0.137 


0.130 


-0.005 


0.001 


0.035 


0.023 


1 


884 


1 


885 


205 +1.63821461(12) 


1 


834 


2 H 


■ +3~ 


-0.236 


0.175 


0.014 


0.022 


-0.003 


-0.048 


1 


785 


1 


791 


203 +1.62225787(12) 


1 


800 


2 H 


" +3" 


-0.363 


0.281 


0.022 


0.027 


-0.003 


-0.063 


1 


737 


1 


748 


201 +1.605(2) 


1 


802 


n- 
2 


" + 3~ 


-0.684 


0.505 


0.049 


0.057 


-0.001 


-0.133 


1 


669 


1 


705 


199 +1.60(2) 


1 


808 


2 H 


■ + 3~ 


-0.649 


0.475 


0.047 


0.057 


-0.003 


-0.130 


1 


678 


1 


712 


197 +1.58(2) 


1 


820 




2+ 


-0.426 


0.332 


0.025 


0.036 


0.0 


-0.070 


1 


750 


1 


764 


195 +1.58(4) 


1 


831 




2+ 


-0.390 


0.303 


0.023 


0.034 


0.0 


-0.065 


1 


766 


1 


777 


193 +1.5912(22) 


1 


838 




2+ 


-0.717 


0.442 


0.076 


0.089 


0.001 


-0.198 


1 


639 


1 


695 


191 1.588(4) 


1 


852 




2+ 


-0.370 


0.281 


0.023 


0.032 


0.0 


-0.065 


1 


787 


1 


798 


189 


1 


905 




2+ 


-0.413 


0.310 


0.027 


0.034 


0.0 


-0.076 


1 


829 


1 


843 


187 


1 


908 




2+ 


-0.391 


0.292 


0.026 


0.031 


0.0 


-0.073 


1 


836 


1 


848 



TABLE IX: PC contributions to the magnetic moments (in fi^) of odd Bi isotopes, Ao = l<7g/ 2 - The notation is explained in 
the text nearby Table IV. 



A 


/^cxp 


Mo 




V 






r^GDD 




) $ Vend 


(5/itot 


Mth 


Ath 


213 


+3.717(13) 


3.659 




2+ 


-0.196 


0.170 


0.013 


-0.013 


-0.001 


-0.014 


3.645 


3.646 


211 


3.5(3) 


3.651 




2+ 


-0.029 


0.026 


0.002 


-0.001 


0.0 


-0.002 


3.650 


3.650 


209 


+4.1106(2) 


3.691 




all 


-0.146 


0.096 


0.033 


0.041 


0.002 


-0.015 


3.619 


3.620 


207 


+4.0915(9) 


3.692 


2 H 


■ + 3~ 


-0.575 


0.485 


0.043 


0.001 


0.002 


-0.045 


3.647 


3.653 


205 


+4.065(7) 


3.696 


2 H 


■ +3~ 


-1.072 


0.922 


0.074 


-0.032 


0.003 


-0.073 


3.624 


3.640 


203 


+4.017(13) 


3.696 


2 H 


■ +3~ 


-2.263 


1.965 


0.151 


-0.109 


0.004 


-0.143 


3.553 


3.607 


201 


4.8(3) 


3.692 


2 H 


■ +3~ 


-2.344 


2.037 


0.158 


-0.109 


0.004 


-0.145 


3.547 


3.604 


199 


+4.6(4) 


3.689 




2+ 


-3.330 


2.925 


0.215 


-0.218 


0.002 


-0.188 


3.501 


3.590 



The calculation scheme is the same as for the Tl and 
Bi chains, in particular, we use the BM values for the 
phonon gyromagnetic ratio. On the base of experience of 
calculations in the lead region, in all superfluid tin nuclei 
we consider the 2f states only except for the magic 100 Sn 
and 132 Sn where we include also the 3j~ state. In Table 
X the characteristics of the 2f states, excitation energies 
and B(E2) values, are given. They, just as those for the 
lead isotopes, were reported in [15|. As it is follows from 
the above analysis, the excitation energy UJ2 is of primary 
importance for accuracy of calculations of the PC correc- 
tions. We see that, with several exceptions, L02 value is 
reproduced with accuracy 100 + 200 keV i.e. sufficiently 
well. The B(E2) values characterize the "strength" of 
the phonon coupling amplitudes g^, and they are also 
reproduced reasonably. 

In Table XI, the Z-factor values are shown for three 
states Ao nearby the Fermi surface, just the same which 
are involved in In and Sb nuclei we analyze. We see that 
there are several cases where the Z-factors are very small, 
1 — Z 1. This is a signal that there is a catastrophic 
term with a small denominator in Eq. PB]) . E.g., in 
the 110 Sn case, for Z(lg7/ 2 ) we have £i g7/2 — £2d 5/2 — 
Wj, = 0.08 MeV. It's clear that in such a situation the 
perturbation expansion in g\, in fact, in (3^ of Eq. (|46[) 
is absolutely non-realistic, and an exact account of the 



TABLE X: Characteristics of the low-lying 2^~-phonons in 
even Sn isotopes, u 2 (MeV) and B(£2,up) x 10 4 (e 2 fm ). 



A 


Ll>2 


exp 


B(E2) th 


B(E2) 


exp 


102 


1 


453 


1 


472 





065 








104 


1 


388 


1 


260 





107 








106 


1 


316 


1 


207 





142 


0.195 


(0 


039) 


108 


1 


231 


1 


206 





155 


0.222 


(0 


019) 


110 


1 


162 


1 


212 





188 


0.220 


(0 


022) 


112 


1 


130 


1 


257 





197 


0.240 


(0 


014) 


114 


1 


156 


1 


300 





193 


0.24 


(0 


05) 


116 


1 


186 


1 


294 





182 


0.209 


(0 


006) 


118 


1 


217 


1 


230 





172 


0.209 


(0 


008) 


120 


1 


240 


1 


171 





152 


0.202 


(0 


004) 


122 


1 


290 


1 


141 





158 


0.192 


(0 


004) 


124 


1 


350 


1 


132 





147 


0.166 


(0 


004) 


126 


1 


405 


1 


141 





120 


0.10 


(0 


03) 


128 


1 


485 


1 


169 





094 


0.073 


(0 


006) 


130 


1 


610 


1 


221 





055 


0.023 


(0 


005) 


132 


4 


327 


4 


041 





104 


0.11 


(0 


03) 


134 


1 


142 





725 





033 


0.029 


(0 


005) 



dangerous terms is necessary. Fortunately, the nuclei for 
which experimental data exist are out of this catastrophic 
region. 
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TABLE XI: Z-factor values of proton levels close to the Fermi 
surface induced by the low-lying phonons in the tin chain. 



A Z(lg 9/2 ) Z(2d 5/2 ) Z(lg 7/2 ) 



100 


0.939 


0.910 


0.933 


102 


0.818 


0.783 


0.524 


104 


0.716 


0.673 


0.344 


106 


0.634 


0.592 


0.238 


108 


0.584 


0.549 


0.158 


110 


0.509 


0.481 


0.103 


112 


0.482 


0.464 


0.458 


114 


0.510 


0.484 


0.238 


116 


0.609 


0.582 


0.414 


118 


0.578 


0.547 


0.454 


120 


0.614 


0.591 


0.607 


122 


0.634 


0.605 


0.591 


124 


0.674 


0.648 


0.651 


126 


0.734 


0.712 


0.725 


128 


0.800 


0.780 


0.798 


130 


0.914 


0.878 


0.836 


132 


0.962 


0.923 


0.968 



Just as for the lead chain above, we display in Fig. 21 
the components of the normal amplitude gf^ in the 118 Sn 
nuclei which is chosen as a representative of the tin chain. 
For graphical convenience the spin-one components are 
multiplied by the factor of 5. As before, we deal with the 
proton vertices only. We see that again this mode is a 
typical surface vibration, the spin-one component is very 
small and can be neglected. For completeness we present 
again in Fig. 22 the neutron spin-zero components of 
the anomalous amplitudes g^) gf\ We see that both 
of them are rather large and have the surface maxima 
of opposite sign. In the result, the amplitude 3™ + = 
+ g^ is very small whereas the amplitude <?^ n '~ = 

g^ — g\^ is comparable with the normal amplitude gf^ ■ 
This point is out of the scope of this article but it worth to 
note that the latter combination determines the collective 
phenomena in superfluid sub-system [l|, H| . 

Let us go to PC corrections to magnetic moments and 
begin with odd In nuclei. In this chain, there are 13 iso- 
topes with known experimental values of magnetic mo- 
ments of the ground state 9/2+ corresponding to the g^J 2 
hole in the corresponding even core of Sn. The theoret- 
ical predictions for the magnetic moments are displayed 
in Fig. 23 whereas separate components of the PC cor- 
rection are given in Table XII. We see that the TFFS 
prediction without PC corrections overestimates the ex- 
perimental value by ~ 0.4 fx^. The value of <5/Xtot, col- 
umn 10 of Table XII, is always negative and in general 
large enough to compensate this excess. In 109 - 113 i n iso- 
topes the PC correction is even too large. This correlates 
with small values of the Z ~ 0.5 which corresponds to 
dYj/de ~ — 1. It is worth noting that even the biggest 
total PC correction 8^i t ot — —0.6 /ijv is the result of an 
almost exact compensation of the two ten times bigger 
corrections Sfi Z and SficGD- Therefore any additional 
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FIG. 21: Components of the normal amplitude g 2 in 118 Sn. 
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FIG. 22: Components of the neutron anomalous amplitudes 
g^ and g™ in 118 Sn. 



term can change the results significantly. The PC cor- 
rection is reduced in this case by almost a factor two if 
we use the ansatz (|54p . fi t h Ath- Remind that it is an 
approximate way to take into account some higher order 
terms in g\ corresponding to partial summation of the 
pole end diagrams. As it was discussed above, so large 
PC effect in the Z-factor is a signal of a dangerous term 
with a small energy denominator in the sum of Eq. (|24[) 
which determines the 5/i z and in its counterpart in the 
sum of Eq. (|34l) for the S^ggd term. More exact ac- 
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FIG. 23: Magnetic moments of In isotopes. 



FIG. 24: Magnetic moments of Sb isotopes. 



count for these dangerous contributions is necessary of 
course. But now we can hope that in so anomalous cases 
the truth is somewhere between /i t h and /ith- 

Let us now go to Sb isotopes. Here 4 experimental 
magnetic moments in the ground states 5/2+ of the odd 
115-121 gk i so topes are known which correspond to the 
2c^ 2 particle states and 7 moments for the lgy 2 state in 
I2i-I33gt) isotopes, in the 121 Sb isotope we dealing with 
the first excited state. The total magnetic moments are 
displayed in Fig. 24 whereas separate PC corrections are 
represented in Table XIII. For the lc% /2 state the TFFS 
without PC corrections describes the data better than 
for In isotopes discussed above. The maximum of the 
disagreement, ~ 0.3^, is for 133 Sb. Unfortunately, the 
PC correction in the vicinity of the magic 132 Sn nucleus, 
just as for neighbors of 208 Pb considered previously, is 
rather small. This correction is rather small also for other 
isotopes, with exception of two left members of the chain. 
Here the term S/i' end is bigger than usual and leads to the 
noticeable value of 6 foot- It is worth to mention that the 
121 Sb nucleus, Ao = I.97/2, is the only case with positive 
value of (5/itot- As usual, the ansatz /i t h — > Mth damps 
the PC effects, in this case resulting in better agreement 
with the data. 

For d 5 /2 states, the disagreement between TFFS and 
the data is significant, about 0.6 /xjv- The size of the PC 
correction (5/itot is of this order or even bigger. Again the 
use of the ansatz (|54p diminishes the PC effect and again 
the experimental magnetic moments are between /i t h and 
/ith values. 

Thus, we found that in semi-magic nuclei the PC cor- 
rections to the magnetic moments are significant. If we 
exclude Bi isotopes where they are small, the PC cor- 



rections have always the sign necessary to improve the 
agreement with the data. 



VIII. CONCLUSION 

Within the self-consistent TFFS, we developed a model 
to calculate the corrections to magnetic moments due to 
coupling to the low-lying phonons in odd magic and semi- 
magic nuclei. The main idea of the model is to consider 
explicitly only those PC diagrams which are sensitive to 
the nucleus under consideration and the state Ao of the 
odd nucleon. The omitted diagrams are supposed to be 
included into the universal TFFS parameters. 

Among semi-magic nuclei we limit ourselves to those 
where the odd nucleon belongs to the non-superfluid sub- 
system. This simplifies the consideration significantly as 
all PC corrections can be calculated without taking into 
account the pairing effects. The perturbation theory up 
to the order of g\ is developed, where is the phonon- 
particle coupling vertex. The method guarantees the to- 
tal angular momentum conservation. For this aim, an 
ansatz is proposed which takes into account the so-called 
tadpole term. 

Three kinds of the PC corrections are considered. The 
end correction, see Fig. 3, is the first one. The main 
(pole) part of the end correction, S/i z , results in \fZ- 
factors on the "ends' of the effective field Vva- The sec- 
ond one denoted as S^ggd, see Fig. 5, describes the 
effect of the induced interaction due to the L-phonon ex- 
change. The third one labeled as 5[Igdd, see Fig. 6, de- 
scribes the effect of the magnetic moment of the phonon. 
This correction contains two terms of different analytical 
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TABLE XII: PC contributions to the magnetic moments (in /ijv) of odd In isotopes, Ao = 159/2, L* = 2 + . The notation is 
explained in the text nearby Table IV. 
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TABLE XIII: PC contributions to the magnetic moments (in /ijv) of odd Sb isotopes. The notation is explained in the text 
nearby Table IV. 
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structure. The term $H GDD i s regular at small phonon 

excitation energy u>l whereas the term S[i GDD is singular 
at ujl — > 0. Each of the three types of PC corrections has 
a tadpole-like counterpart and the one shown in Fig. 7, 
corresponding to the S/igdd term, is obviously the most 
important one. In particular, it possesses the same sin- 
gularity as the 5^q DD term. We suppose that the sum of 

tadpole-like terms and the <5Mgdd one compensate each 
other. It is motivated by the observation that such an 
ansatz leads to the total angular momentum conservation 
with g\ accuracy. 

In magic nuclei, the main PC correction comes from 
the 3^ phonon, whereas in non-magic ones the 2 + - 
phonon dominates. The main observation of our cal- 
culations is the very strong cancelation of the first two 
corrections, Sfi z and <5mggd, The scale of the PC cor- 
rections may be characterized with the quantity (1 — Z) 
which determines directly the term 8fi z . With only one 
exception, the total correction to the magnetic moment 
(J^tot has the same sign as 5[i z , and the absolute value 
is up to ten times smaller than the individual contribu- 
tions. As the result of this cancelation, in magic nuclei 



where the value of (1 — Z) ~ 0.1 is typical, PC correc- 
tions are as a rule negligible. In any case, it is so for 
all magnetic moments where experimental data exist. In 
particular, PC corrections can not help to solve the old 
for the TFFS theory problem of the magnetic moment of 
the ground state of 209 Bi where the theoretical value de- 
viates by ~ 0.5 (in from the experiment. Evidently, the 
standard form of the spin-dependent LM amplitude is 
not sufficient in this case and additional terms including 
the spin-orbit and tensor in the cross-channel interaction 
should be taken into account. Thus the problem should 
be solved at the RPA level. 

Simultaneously we calculated the gyromagnetic ratios 
<?£ h of all low-lying phonons in 208 Pb. For the 3j~ state 
it is rather close to the BM model prediction g£ h = Z/A 
which is a result of its strong collectivity. Other L- 
phonons are much less collective and their gyromagnetic 
ratios strongly differ from the BM ones. 

We calculated the PC corrections to four chains 
proton-odd semi-magic nuclei. The odd Tl and Bi iso- 
topes are the odd neighbors of the even lead isotopes. 
The Bi chain contains 7 nuclei with known magnetic mo- 
ment in the same l/ig/ 2 ground state and the PC cor- 
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rection Sfitot is much less than necessary to explain the 
data in non-magic isotopes too. Evidently, the problem 
of 209 Bi should be solved first at the RPA level which 
will reduce the disagreement for the whole Bi chain. For 
Tl isotopes there are 10 isotopes with known magnetic 
moment in the 3s7^ ground state. For the magic isotope 

207 T1 TFFS agrees with the data perfectly well and the 
PC correction as it was said above is negligible. For non- 
magic Tl isotopes the typical disagreement is ~ 0.2 //at 
and PC corrections reduces it by one half. For odd- 
proton neighbors of the even tin isotopes, the odd In 
and Sb chains, the PC effects are stronger than in the 
lead chain. In such a situation, sometimes the use of 
the perturbation theory in g 2 L seems questionable and we 
suggested the ansatz given in Eq. (|54|) of an approximate 
account for the higher order terms in g\. It corresponds 
to partial summation of the pole end diagrams in Fig. 3 



itself and in the "ends" of Fig. 5 and Fig. 6 as well. This 
ansatz, /i t h — > fhb, reduces the PC effect that is impor- 
tant in several cases. In general, in both the chains the 
PC corrections improves the agreement. 
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